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In this article I will review the theory behind the gravitational wave driven r-mode 
instability in rapidly rotating neutron stars and discuss which constraints can be derived 
from observations of spins and temperatures in Low Mass X-ray Binaries. I will discuss 
how a standard, ‘minimal’ neutron star model is not consistent with the data, and 
discuss some of the additional physical mechanisms that could reconcile theory with 
observations. In particular I will focus on additional forms of damping due to exotic cores 
and on strong mutual friction due to superfiuid vortices cutting through superconducting 
flux tubes, and examine the repercussions these effects could have on the saturation 
amplitude of the mode. Finally I will also discuss the possibility that oscillations due 
to r-modes may have been recently observed in the X-ray light curves of two Low Mass 
X-ray Binaries. 
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1. Introduction 

Rapidly rotating Neutron Stars (NSs) are one of the main targets for Gravita¬ 
tional Wave (GW) detectors such as Advanced LIGO and Virgo,® as there are 
several mechanisms that could lead to non-axisymmetric deformations of the star 
and ultimately to a continuous GW signal.® In particular there are several modes 
of oscillation of the NS that can result in GW emission, and in this context the 
so-called r-mode has attracted considerable attention, as it is generically unstable 
to GW emission,®® and could thus grow to large amplitudes and offer the best 
detection prospects. 

The possibility of detecting GWs from NS modes of oscillation is particularly 
exciting as gravitational wave asteroseismology would allow us to probe the interior 
structure and composition of NSs in great detail, in the same way as electromag¬ 
netic asteroseismology has, in recent years, allowed us to significantly enhance our 
understanding of many other stars, from white dwarfs to red giants.®® 

Let us remind the reader that, with a mass roughly equal to that of the sun 
compressed in a 10 km radius, NSs are one of the most compact objects in the 
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universe and their interior density can easily surpass nuclear saturation density. 
Furthermore NSs are cold objects as, although their core temperatures are in the 
region of T ss 10® K, the thermal energy is small compared to the Fermi energy 
in most of the star. This can have important consequences for the dynamics of 
the system, as neutrons will pair and form a superfluid, while protons will be su¬ 
perconducting. Probing the interior structure of these objects would allow us to 
probe physics in an entirely different regime from ground based experiments, such 
as heavy ion colliders, which generally probe the low density and high temperature 
section of the QCD phase diagram.l^ While at asymptotically high densities quarks 
are though to pair in the so-called Colour-Flavour-Locked (CFL) phase,!^ for real¬ 
istic NS densities the ground state of matter is unknown, and only astrophysical 
observations will be able to solve this problem. 

In order to obtain astrophysical constraints and set the theoretical basis for 
the analysis of future GW data, in this article we will focus on some of the most 
promising systems for GW detection: Low Mass X-ray Binaries (LMXBs). These 
are binary systems in which a compact object (in our case a neutron star) is ac¬ 
creting matter from a less evolved companion that fills its Roche lobe. As matter is 
transferred from the secondary to the primary it forms an accretion disc and even¬ 
tually is accreted by the NSs, spinning it up. This is thought to be the mechanism 
by which old, long period, stars are recycled to millisecond periods and millisec¬ 
ond radio pulsars are eventually formed.l^ The reason LMXBs are invoked as GW 
sources lies in an observational puzzle: while accretion should be able to spin-up the 
NS to its Keplerian break-up frequency (which is equation of state dependent, but 
generally above 1.5 KHz), there appears to be a cutoff at around 700 Hz in the dis¬ 
tribution of observed spins for both the LMXBs and the millisecond radio pulsars. 
It was thus suggested that GW emission could provide an additional spin-down 
torque that would balance the torque due to accretion and stall the spin-up.l^^lii^ 
Several mechanisms have been proposed, including ‘mountains’ supported by the 
crust, ^ the cor^i^E^or confined by the magnetic field,E^and unstable modes of os- 
cillation.l^lS Although detection will be challenging for many of these scenarioJi^^ 
and in several cases it is more likely that the disc/magnetosphere interaction is 
dominating the torques these systems still allow for the best constraints on 
the physics of the r-mode instability and the interior dynamics of NSs. Further¬ 
more recent detections of oscillations in the X-ray light curve of two LMXBs may 
be interpreted as r-modes perturbing the electromagnetic emission, and open the 
fascinating prospect of combining GW and electromagnetic signals to study NS 
interiors 

2. The r-mode instability window 

We begin by focusing on the r-mode instability. An r-mode is a fluid mode of 
oscillation for which the restoring force is the Goriolis force. It thus only exists in a 
rotating star. In Newtonian gravity, and to first order in the rotational frequency 
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of the star, it is purely toroidal and for the Eulerian velocity perturbation Sv one 
has (for an in depth review and a discussion of relativistic effects see e.g. references 
[2T] and [22]): 


Sv = a 



mY 


B „icot 
Im^ 


( 1 ) 


where = [l{l + l)]-^/^rY x (rVF/™) is the magnetic-type vector spherical 

harmonic (with Yim the standard spherical harmonics), R is the stellar radius and 
a the dimensionless amplitude of the mode,^ while the frequency of the mode uj 
takes the form!^ 

‘^TYl 


This mode is interesting for our discussion because it is generically unstable to GW 
emission. In particular if we examine the pattern speed for an r-mode this is, in a 
frame rotating with the star: 


(7 •y' - 


UJr 

m 


212 


( 3 ) 


In the inertial frame on the other hand one has 


CTi = 


a-i)(^+2) 
/(/ + 1 ) 


12 


( 4 ) 


so that a mode that is retrograde in the rotating frame appears prograde in the in¬ 
ertial frame. The r-modes thus satisfy the criterion for the so-called Chandrasekhar- 
Friedman-Schutz (CFL) instability,l^l^ which allows for the star to find lower en¬ 
ergy and angular momentum configurations in which the mode amplitude can grow. 
Note that other modes of oscillation can also be unstable, and another candidate 
for GW detections is the f-mode.^^ This mode is not, however, generically unstable, 
but only goes unstable above a critical frequency. Furthermore in cold systems, 
such as the LMXBs we consider, the f-mode instability is generally stabilised by 
viscosity due to superfluid mutual friction.l^^H^ The strongest contribution to GW 
emission is due to the I = m = 2 r-mode an in this case for an n = 1 polytrope 
(which we shall take as our equation of state for all the following estimates) the 
growth time of the instability is!^ 


Tgw = -47 
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Naturally the mode amplitude can only grow provided viscosity cannot damp the 
instability faster than GW emission can drive it. Viscosity in neutron star interiors 
is not entirely understood, and it is thus difficult to model in detail. To illustrate the 
problem we will thus first introduce a ‘minimal’ NS model, corresponding the most 
commonly considered setup, and calculate the region of parameter space in which 
the r-mode can grow. In the following we will then show that X-ray observations 
of spins and temperatures of NSs in LMXBs point towards the need for additional 
physics in our model. 
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In our ‘minimal’ model we will assume that, at high temperatures, bulk vis¬ 
cosity due to modified URCA reactions provides the main damping mechanism, 
while at low temperatures the main contribution is from shear viscosity, due to 
standard scattering processes (mainly electron-electron in superfluid mattei!^, or 
from viscosity at the crust-core interface. We assume that there are no dynamically 
important magnetic fields or superfluid degrees of freedom, no exotica in the core, 
and, as previously stated, we take an n = I polytrope as equation of state. With 
these assumptions we can calculate the damping timescale associated with bulk 
viscosity, , and shear viscosity, Tgv 


Thy = 2.7 X 10 
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10 km 

where M is the mass of the star, R is the radius, P is the rotation period and T 
is the core temperature. For damping due to Ekman pumping at the crust/core 
interface we use the estimate of Glampedakis & Anderssor(22ll2U with a slippage 
parameter S = 0.05. The slippage parameter is essentially the ratio between the 
curst/core velocity difference and the mode velocity, and accounts for the fact that 
the crust is not completely rigid, but can participate in the oscillation. In this case 
the damping timescale is: 


Tefc = 3 X 10“^ 
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In order to understand how a system evolves in the presence of an unstable 
r-mode it is useful to consider an ’instability’ window, i.e. a region in parameter 
space in which the mode is unstable. In order to do this we will fix the mass of 
the star at M = I.4M0 and the radius at i? = 10 km and examine the region in 
the spin frequency vs temperature plane in which the mode can go unstable. The 
boundary of this region corresponds to the points in which the damping and driving 
timescales are equal, i.e. to the solutions of 


I 


Tgw 



(9) 


where Tyi is the viscous damping timescale for process i acting in the star. The 
result is shown in figure [TJ 

Phenomenological evolution equations for the frequency of the star, the mode 
amplitude a and the thermal energy Ex of the star take the formP^ 
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dt Tv 1 + a^Q 

da a a 1 — a^Q 

dt Tgyy Ty I P Q 

dEx 


( 10 ) 

( 11 ) 


dt 


El, + Ey 


( 12 ) 






















September 16, 2015 0:22 WSPC/INSTRUCTION FILE haskellrmode 


R-modes in neutron stars 5 


where ry is the damping timescale due to the dominant viscous mechanisms, the 
dimensionless parameter Q = 9.4 x 10“^ for an n = 1 polytrop^^ and E^, is the 
luminosity due to the modified URCA process!^ 


= 1.1 X 10' 
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with p the density, while Ey is the reheating due to shear viscosity as the mode 
growspil 


— = 3.3 X 10-2- 

dt Tgy 


(14) 


A typical accreting NS, with a core temperature of around 10® K will be spun 
up by accretion into the unstable region. The r-mode then rapidly grows to large 
amplitude, resulting in fast heating. Eventually the thermal runaway is halted by 
cooling due to neutrino emission and the star simply spins down due to GW emis- 
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Fig. 1. The r-mode instability window for a 1.4 Mq, R = 10 km NS. The equation of state is an 
n = 1 polytrope and we consider the ‘minimal’ model described in the text. In the shaded region 
the mode is stable, and unstable above it. At low temperatures the main source of damping is 
given by Ekman pumping the crust-core boundary, while at high temperature bulk viscosity gives 
the main contribution. We also show, both for small and large amplitudes a of the mode, the cycle 
that an accreting system would follow in the temperature-frequency plane. 


However the above equations are only valid as the mode is growing, and in 
reality the amplitude will saturate at a value as << 1 due to non-linear couplings 
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to other modes. We can see from equation [14] that this saturation amplitude plays 
a critical role in determining the amount of heating and ultimately how far into the 
instability window a NS can move. 

If the saturation amplitude is very large (og ~ I) the system will move well 
into the unstable region, but the evolution will be fast and the duty cycle very 
short,less than Ri 1% . Conversely, if the mode saturates at low amplitude 
(a « 10“®), as suggested by calculations of non-linear couplings to other modes, 
the duty cycle is much longer but the system will remain close to the instability 
curve.l^ In both scenarios it is highly unlikely to observe a system far inside the 
unstable region. 

3. Observati 



Fig. 2. Comparison between the observed spins and core temperatures of NS and the r-mode 
instability window in our ‘minimal’ model, for a M = 1.4Mq and i? = 10 km star, described by 
an n = 1 polytrope. The error bars are due to the uncertainty on the composition of the outer 
layers of the star. There is clearly a large number of sources in the unstable region, while our 
expectation is that they would rapidly be spun out for large saturation amplitudes of the mode, 
or remain close to the instability curve if as is small. 


The main conclusion of the previous section is that one would not expect to 
see systems in the instability window, and that they should either be close to the 
instability curve or well below it. We can test this prediction for our ‘minimal’ 
model by populating the window with observations of spins and temperature in 
LMXBs. The spins of NSs in LMXBs can either be measured directly in systems 
for which coherent X-ray pulsations are detected, or inferred from the frequency of 
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oscillations seen in thermonuclear type I X-ray burstsP^^^ Temperatures, on the 
other hand, require some degree of modelling, given that we are interested in the 
core temperature. Surface temperatures of NSs can be estimated from black body 
fits to spectra of LMXBs in quiescence. To obtain the core temperature one can 
then assume that the interior is roughly isothermal, and model the exterior layers of 
the NS to obtain a relation between the surface temperature and the temperature 
at the base of the envelope.^ 

In figure [2] we see the result of populating our ‘minimal’ instability window for 
a M = I. 4 M 0 and i? = 10 km NS, with data from LMXBs, as obtained by Haskell, 
Degenaar and Ho.^^ The error bars in the figure are due to the uncertainty in the 
composition of the outer layers. It is quite obvious that there is a large number 
of sources inside the instability window, in a region that would not be permitted. 
Given the limited sample of systems for which one has measurements of both spins 
and temperatures (22 in the case of figure [5]) one would not expect to find any in 
the unstable region.l^ The conclusion is quite robust and was found to hold even for 
different masses and equations of state, 1^1 and to be consistent with more detailed 
modelling of the core temperatures of LMXBs.l^ 

This problem clearly shows that the ‘minimal’ model we described, and that 
is often considered in NS physics, needs to be re-evaluated and additional physics 
must be included. In the following we describe some of the main mechanisms that 
could make our theoretical understanding consistent with observations. 


4. Additional damping mechanisms 

First of all we examine the possibility that additional physics, such as exotic par¬ 
ticles in the core or strong superfluid mutual friction, may provide strong damping 
at low temperatures, and modify the shape of the instability window. Let us thus 
outline some of the most promising mechanisms. 


4.1. Mutual friction 

Neutrons in the interiors of NS are thought to pair and form a large scale super¬ 
fluid condensate that can oscillate independently from the proton-electron fluid, to 
which it is only very weakly coupled.!^ R-modes in superfluid neutron stars have 
been studied in detail by several authors,l2^^^ who have found that two families of 
modes can now exist, one in which the fluids are mainly co-moving and another 
in which neutrons and protons are counter moving which, however, only exists in 
non-stratified stars.^^^ Rotation ‘mixes’ the modes and, to second order in rotation, 
the standard co-moving r-mode also has a counter-moving component which, for 
millisecond spin-periods such as those we are considering, can be quite large and 
lead to damping via vortex mediated superfluid mutual friction. The strength of the 
mutual friction is usually quantified by a dimensionless parameter TZ which encodes 
the microphysics that gives rise to the effect. The coupling timescale between the 
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superfluid neutrons and the electrons then scales as r ~ l/2fllZ for IZ « 1 as is 
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Fig. 3. The r-mode instability window for strong mutual friction. There is a considerable difference 
in the windows depending on the superfiuid pairing gaps that are used: we show the ’strong’ (left) 
and ’weak’ (right) models described in Haskell et al. 2009.1^ In both cases the theoretical model 
agrees with observations for TZ ~ 0.01, which is in the possible range for the vortex flux tube 
cutting mechanism. 


The most commonly considered mutual friction mechanism is the scattering 
of electrons of magnetised vortex cores,which gives TZ ~ 10“^. In this case 
the stabilising effect on the r-mode is weak and the instability window is not al- 
tered.l^^ It is thus interesting to ask how large the parameter TZ needs to be for 
the instability window to be consistent with observations. As can be seen in fig¬ 
ure [3] the instability curve depends quite strongly on the superfluid pairing gaps 
that are used (which can also have an impact on the strength of shear viscositjl^, 
but is still generally consistent with observations for TZ w 10“^. This value could 
be consistent with what is expected if the protons in the outer core form a type 
II superconductor and the superfluid vortices ‘cut’ through superconducting flux 
tubes, ^ in which case one expectsP^ 
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where B is the strength of the magnetic field. The exterior dipole component of the 
magnetic held in LMXBs is generally inferred to be in the region B rs 10® — 10® 
G, so the interior held would have to be signihcantly stronger for this mechanism 
to work. It is important to note though that in equation (I15|) the mutual friction 
coefficient TZ is velocity dependent, and we shall see in the following that this can 
have important consequences for the saturation amplitude of the r-mode. 
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4.2. Exotica in the core 

If exotic particles, such as hyperons or deconfined quarks, are produced in the core 
of the NS this generally leads to an increase of the bulk viscosity at low temperatures 
which significantly alters the instability window. In figure |4] we show the instability 
window in the case of a strange star with M = IAMq and i? = 10 km.l^ The 
main point to note is that the increase in bulk viscosity at around T Ri 10® K 
can explain the presence of several systems in this region. The window depends 
slightly on the model that has been chosen for the shear and bulk viscosity, and 
on the parameters of the theory, however Alford and collaborator^^ have shown 
that the main features of the instability window are remarkably insensitive to the 
exact microphysical description of matter, and models with a sizeable ungapped 
quark are generally consistent with observations.^^ Note, however, that if quarks 
in the core are paired in the CFL phase, viscosity is much lower and the instability 
window is essentially that of the ‘minimal’ model.^ The appearance of hyperons 
in the core can also lead to an increase in bulk viscosity at low temperatures.^^l^ 
If r-modes in LMXBs are stabilised by strong bulk viscosity at low temperatures, 
this poses, however, a new theoretical challenge, as one would expect the systems 
to cool and spin down along the instability curve, after accretion ceases. This is 
at odds with the observations of several fast (with spin frequencies up to 700 Hz) 
millisecond radio pulsars. This issue is not present if shear viscosity is enhanced at 
high temperatures, which may be the case in the presence of pion condensation.^^ 


4.3. The crust/core interface 

In our ‘minimal’ model we have included the effect of damping due to a viscous 
boundary layer at the crust/core interface. We have accounted for the possibility 
that the crust may not be completely rigid but participate in the oscillation by 
introducing a slippage parameteJ^I^ S that is essentially the ratio between the 
jump in velocity at the crust core interface and the mode velocity. For small values 
of S the crust can oscillate freely, while for 5 = 1 the crust is completely rigid, 
and in this case all the observed systems would be r-mode stable. In our model 
we have set S = 0.05, however the crust does not respond in the same way at 
all frequencies, as resonances between the r-mode and torsional oscillations of the 
crust are possible. This leads to the possibility that the rigidity of the crust, i.e. 
the parameter S, may be frequency dependent,S^l^and more detailed modelling of 
the coupling between crust and core is necessary to understand if this mechanism 
is consistent with observations of temperatures and spins in LMXBs. 

Note, however, that if the crust/core transition is not sharp, but more gradual 
and proceeds via several phase transitions that give rise to the so-called pasta 
phases, then viscosity would be much weaker and by mainly due to standard shear 
viscosity .1^ 
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Fig. 4. The r-mode instability window in the case of a strange quark star with M = 1.4M0 and 
it = 10 km.EniThe exact shape of the window depends slightly on the choice of model parameters, 
such as the value of the strong coupling constant Os and the mass of the strange quark rris but 
the qualitative features, namely the increase of viscosity around T = 10® K are fairly insensitive 
to the exact microphysical description of the quark core.lSl 


4.4. Interactions with superfluid modes 

Another interesting possibility that was introduced by Gusakov and collaborator^^ 
is that of interactions between ‘superfluid’ inertial modes and the standard r-mode, 
which can take place at fixed ‘resonance’ temperatures and lead to ‘spikes’ in the 
instability window as the mode changes characterThis is due to the fact that 
the ‘superfluid’ inertial modes are counter-moving oscillations of the fluid already 
to leading order in the rotation rate of the star il. For these ‘superfluid’ modes 
one has for the counter-moving component of the velocity Wg « 0(fl), compared 
to w„ Ri for the ‘normal’ r-mode. As counter-moving motion is the main 

driver behind mutual friction, a large w at leading order in rotation leads to a 
short mutual friction damping timescale, such that for a ‘supefluid’ mode one has 
T^p ~ 0(n), rather than r^jp ~ 0(12)^ for a normal r-mode. Mutual friction is 
the main viscous process acting on ‘superfluid’ inertial modes and provides efficient 
damping, while the coupling to GW emission is much weaker.l^The r-mode is thus 
stabilised by the interaction with these modes. 

This scenario makes an observational prediction, namely that one should un¬ 
cover a population of HOt and Fast Non-Accreting Rotators (HOFNARs) that are 
essentially ‘hot windows’, i.e. systems in which the NSs has been re-heated by the 
r-mode instability and appears as a thermally emitting isolated system after the 
accretion phase is over. ^Future X-ray observations could thus confirm the viability 
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of this scenario. 

5. Saturation amplitudes 

Up to now we have described mechanisms that can enhance the viscosity in the 
temperature and frequency range of LMXB observations, and modify the instability 
window so that these systems are r-mode stable. The ‘minimal’ model for this 
instability window could, however, be made consistent with observations if the 
saturation amplitude of the r-mode is so small that the instability is indeed present 
in all the observed systems, but at such a low-level as to not impact on the spin or 
thermal evolution of the stars. 

This scenario requires very small amplitudes as ~ 10“® — lO“®,SnillI] generally 
much smaller than the values obtained for saturation due to non-linear couplings 
with other inertial modes,!^ which lead to saturation amplitudes Os ~ 10“®. It has 
been suggested^ that if the crust/core transition is smeared out, for example due to 
the appearance of pasta phases, as previously discussed, then viscosity will be much 
weaker and allow for inertial modes to grow and saturate the r-mode at amplitudes 
tts ^ 10“®. These calculations ignore, however, the effect of mutual friction which 
can quite efficiently halt the growth of inertial modes coupled to the r-mode. 

If the core of the star is in a type II superconducting state, however, mutual 
friction due to superfluid vortices cutting though superconducting flux tube may 
effectively saturate the r-mode at low amplitudes.^ As we have already discussed 
if vortices can cut through flux tubes this leads to strong mutual friction, with TZ Ri 
2.5 X this, however, is not always the case, as the large energy penalty 

of cutting leads to vortices ‘pinning’ to flux tubes for counter-moving velocities Wp 
less than 

/ n \ 1/2 

|w|p « 1.5 X j 

where B is the macroscopic core magnetic field. If vortices are pinned and cannot cut 
through flux tubes there is essentially no mutual friction damping and the r-mode 
can grow (provided other sources of viscosity are not suppressing the instability). 
As the mode grows, however, the amplitude of the counter-moving component of 
the velocity w also grows, according to 

|w|«»A,»(l)'(|-) RU (17) 

with rix the Keplerian breakup frequency of the star and Aq a spin independent 
parameter that lies in the range Aq « 0.1 — When a grows large enough it is 

thus possible to enter the regime in which w > Wp and vortices are forces through 
flux tubes. At this point strong mutual friction due to cutting will take over and 
damp the mode, effectively halting the growth of the instability at an amplitude 
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with ly the spin frequency of the star in Hz. In many cases this amplitude can 
be smaller than the limit set by non-linear mode couplings.^ Even lower satura¬ 
tion amplitudes of a ~ 10“^® are possible in hybrid stars due to periodic phase 
conversion at the interface with the exotic core.l^ 

Another possibility is that the r-mode oscillation will wind up the interior mag¬ 
netic field of the star to the produce a strong toroidal component, which could then 
rapidly suppress the instability if the internal magnetic field of the star is of order 
B I0 ^° gUMII] Note that this mechanisms is expected to be active only if the 
r-mode is driven unstable by GW emission,!^!^ and if the amplitude is not growing 
the effect of the magnetic field will be negligible for the field strengths expected in 
LMXBs.l^SlZll 

The small saturation amplitudes predicted in this section would all lead to GW 
emission well below the detectability level for current and next generation GW 
detectors, for which only emission at the level required to balance the accretion 
torque and explain the observed spin periods would be possible, although chal¬ 
lenging, to detect.li^ Emission at that level would, however, be inconsistent with 
the observed temperatures of LMXBs, as it would heat the stars up more than is 
observed.SSlSal It is thus likely that GWs from r-modes in accreting neutron stars 
in LMXBs will be very challenging to detect, with young NSs being a much more 
promising target.^^ 


6. Electromagnetic observations of r-modes 

Recently Strohmayer & Mahmoodifai^l^ have conducted targeted searches for 
oscillations in the X-ray light curves of several LMXBs and discovered candidate 
features in two of them, XTE J1751-305 and 4U 1636-536. In table [1] we show the 
oscillation frequencies and characteristics of the two systems. 

For the first, XTE J1751-305 the oscillation frequency has been interpreted as 
the rotating frame frequency of the mode, as it is thought that we are observing 
modulations of the X-ray emitting hotspots due to a mode on the surface of the 
star.l^ The best candidate for this mechanism is a surface g-mode,l^^ as it can 
produce large modulations. However r-modes are also likely candidates, especially 
due to the fact that in a stratified star there are interactions (avoided crossings) 
between the g-modes and inertial modes such as the r-modes, that can modify the 
eigenfunctions and the nature of the oscillations.^ If the oscillation is interpreted 
as a global r-mode the amplitude required to explain the observed modulation of 
the light curve is a ~ 10“^ which is very large, and would lead to the system spin¬ 
ning down due to GW emission, which is not observed.^ This problem is somewhat 
alleviated if one accounts for a rigid crust, in which case interactions with crustal 
oscillations can amplify the amplitude at the surface by up to two orders of magni¬ 
tude leading to amplitudes in the bulk of the star of a ~ 10“®. This amplitude 
would still be too large to be consistent with the observed temperature of the star,Sil 
which would require a ~ 10“®, but suggests that further detailed modelling is nec- 
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Table 1. The spin frequency of the star (i/) and observed mode 
frequency for the two systems in whic h an o scillations was de¬ 
tected by Strohmayer &; MahmoodifarP^I^ We also show the 
ration k between the mode frequency and the spin frequency. 


System 

u (Hz) 

Mode frequency (Hz) 

Ratio K 

XTE J1751-305 

435 

249.33 

0.5727597 

4U 1636-536 

582 

835.6440 

1.43546 


essary. 

In the second system, 4U 1636-536, the oscillation appears during a super burst, 
which is thought to be a thermonuclear burst due to unstable burning of a carbon 
layer formed by the ashes of a regular type-I X-ray burst. In this case, consistently 
with the idea that the burning layer spreads rapidly and the emission is coming 
from the whole surface, the observed frequency is interpreted as the inertial frame 
frequency of the mode. As for the previous system, several modes could oscillate at 
the observed frequency, with g-modes and r-modes being the main candidates. 4U 
1636-536, however, is not observed as an X-ray pulsar, and there is no measurement 
for its long term spin-down rate. One cannot thus rule out the possibility of a large 
amplitude r-mode spinning down the star. 

Finally note that both these events took place before the current LIGO science 
runs (in 2002 and 2001 respectively). If such an event were to repeat itself while 
advanced detectors are taking data a large amplitude r-mode would be detected.l^ 
Needless to say a simultaneous X-ray and GW observation of a NS oscillation would 
allow us an unprecedented insight into the interior dynamics of these systems and 
inaugurate the era of NS asteroseismology. 

7. Conclusions 

In this article I have reviewed the main theoretical aspects of the GW driven r- 
mode instability in rapidly rotating neutron stars and described what the instability 
window (i.e. the region in the frequency vs. core temperature plane in which the 
r-mode is unstable) is predicted to be in a ‘minimal’ NS model that assumes a core 
of neutrons, protons and electrons, with no exotica (such as deconfined quarks), 
dynamically significant superfluid degrees of freedom or magnetic fields. 

The predictions of such a model can be compared to observations of spins and 
temperatures in LMXBs. The expectation is that no system should fall inside the 
unstable region. If the saturation amplitude of the mode is small (a m 10“®) as 
predicted by calculations of nonlinear couplings to other modes,^ then the NS 
will never depart significantly from the instability curve. On the other hand, if the 
saturation amplitude is large (a « 1) the system can enter well into the instability 
region, but the timescale to exit it will be very short, leading to a very low duty 
cycle and probability of observing a system in this phase. 

In section|3]I show that the predictions of the ‘minimal’ model are not consistent 
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with observationsas many observed systems would fall in the unstable region. 
Despite uncertainties on the mass of the NS and on the composition of the outer 
layers, this qualitative conclusion is robust and indicates the need to include addi¬ 
tional physics in our model. There are essentially two theoretical solutions to this 
observational puzzle: either there are additional physical mechanisms giving rise to 
strong viscosity at low temperature, and making the observed systems stable, or 
the systems are indeed unstable, but the r-mode saturates at such a low amplitude 
that it does not impact on the spin or thermal evolution of the system. 

Additional viscosity could be due to hyperonJ^^^ or deconfined quarks in the 
core ,S2lI^ strong mutual friction due to vortex/flux tube cutting,^^ resonances with 
crustal torsional mode^^^^ or to interactions between the r-mode and superfluid 
inertial modes.^^ Very small saturation amplitudes, on the other hand, could be 
possible if the growth of the mode is halted before non-linear couplings set in, by 
processes such as vortex/fLux tube cutting,!^ periodic phase transitions in hybrid 
star^^ or by winding up a strong toroidal component of the magnetic field.^ 

Finally I have reviewed recent observations of two oscillation modes detected in 
the X-ray light curve of the LMXBs XTE J1751-305 and 4U 1636 - 536 .^^ 21120 ] both 
cases several modes could lead to the observed frequencies, with the most likely be¬ 
ing surface g-modes or r-modes. For the first system the observed frequency is likely 
to correspond to the rotating frame frequency of the mode, as the surface oscilla¬ 
tion perturbs the emitting hotspot. For the second system, however, the oscillation 
is detected during a thermonuclear super burst, in which burning is likely to be 
occurring over the whole surface of the star. In this case the observed frequency 
would correspond to the inertial frame frequency of the mode. 

Both modes can be interpreted as r-modes, however to explain the observed 
modulation of the X-ray fluxes the corresponding amplitudes would have to be large 
(a ~ 10“^). Although the surface amplitude can be amplified by up to a factor of 
100 by interactions with crustal modes,l^the amplitude is still unrealistically large 
for XTE J1751-305, as the star would have to be hotter than observed, and spin- 
down faster than observed.^In the case of 4U 1636-536, however, a large amplitude 
r-mode cannot be excluded. 

If such an event were to repeat itself while advanced GW detectors are operating, 
a large amplitude r-mode could be detected, allowing for the tantalising possibility 
of a coincident electromagnetic and GW detection, and truly opening the era of NS 
asteroseismology. 
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